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ABSTRACT 

The final stage of terrestrial planet formation is known as the giant impact 
stage where protoplanets collide with one another to form planets. So far this 
stage has been mainly investigated by iV-body simulations with an assumption of 
perfect accretion in which all collisions lead to accretion. However, this assump- 
tion breaks for collisions with high velocity and/or a large impact parameter. 
We derive an accretion condition for protoplanet collisions in terms of impact 
velocity and angle and masses of colliding bodies, from the results of numerical 
collision experiments. For the first time, we adopt this realistic accretion condi- 
tion in iV-body simulations of terrestrial planet formation from protoplanets. We 
compare the results with those with perfect accretion and show how the accretion 
condition affects terrestrial planet formation. We find that in the realistic accre- 
tion model, about half of collisions do not lead to accretion. However, the final 
number, mass, orbital elements, and even growth timescale of planets are barely 
affected by the accretion condition. For the standard protoplanetary disk model, 
typically two Earth-sized planets form in the terrestrial planet region over about 
10 8 years in both realistic and perfect accretion models. We also find that for 
the realistic accretion model, the spin angular velocity is about 30% smaller than 
that for the perfect accretion model that is as large as the critical spin angular 
velocity for rotational instability. The spin angular velocity and obliquity obey 
Gaussian and isotropic distributions, respectively, independently of the accretion 
condition. 

Subject headings: planets and satellites: formation — methods: numerical 
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INTRODUCTION 



It is generally accepted that the final stage of terrestrial planet formation is the giant 
impact stage where protoplanets or pla netary embryos f ormed by oligarchic g rowth collide 



with one another to form planets (e.g. 



Wetherill 



1985 



Kokubo fc Ida 



19981 ). This stage 



has been mainly studied by iV-body simulations. So far all iV-body s imulations hav e 



assumed perfect accretion in which all collisions lead to accretion (e.g., 



Chambers 



Agnor et al. 



1999; 



20011 ). However, this assumption would be inappropriate for grazing impacts 



that may result in escape of an impactor or 
Hydrodynamic (SPH) collision simulations, 



iit-and-run. By per f ormin g Smoothed-Particle 



Agnor fc Asphaug 



(120041 ) estimated that 



more than half of all collisions between like-sized protoplanets do not simply result in 
accumulation of a larger protoplanet, and this inefficiency lengthens the timescale of planet 
formation by a factor of 2 or more, relative to the perfect accretio n case. The accretion 



inefficiency can also change planetary spin. 



Kokubo fc Ida 



(120071 ) found that under the 



assumption of perfect accretion, the typical spin angular velocity of planets is as large as 
the critical spin angular velocity for rotational instability. However, in reality, the grazing 
collisions that have high angular momentum are likely to result in a hit-and-run, while 
nearly head-on collisions that have small angular momentum lead to accretion. In other 
words, small angular momentum collisions are selective in accretion. Thus, the accretion 
inefficiency may lead to slower planetary spin, compared with the perfect accretion case. 

The goal of this paper is to clarify the statistical properties of terrestrial planets formed 
by giant impacts among protoplanets under a realistic accretion condition. We derive an 
accretion condition for protoplanet collisions in terms of collision parameters, masses of 
colliding protoplanets and impact velocity and angle, by performing collision experiments 
with an SPH method. We implement the realistic accretion condition in A^-body simulations 
and probe its effect to further generalize the model of terrestrial planet formation. We 
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derive the statistical dynamical properties of terrestrial plane ts from results o 
of iV-body simulation s and compare the results with those in 



Kokubo et al. 



Kokubo fc Ida 



a number 



2006|) and 



(120071 ) where perfect accretion is adopted. 



In section 2, we outline the initial conditions of protoplanets and the realistic 
accretion condition. Section 3 presents our results, where we show the statistics of collision 
parameters and basic dynamical properties of planets. Section 4 is devoted to a summary 
and discussions. 



METHOD OF CALCULATION 



2.1. Initial Conditions 



We perform iV-body simulations of terrestrial planet formation starting from 
protoplanets. We consider gas-free cases without giant planets residing outside the 
terrestrial planet region to clarify the basic dynamics. To compare with the results with 
perf ect accretion, we adop t the same protoplanet system as those in 



Kokubo et al. 



and 



Kokubo fc Ida 



fl2006f ) 



(120071 ). which is formed by oligarchic growth from a planetesimal disk 



whose surface density distribution is given by 



10 



1AU 



-3/2 



-2 



with inner and outer edges, r- m = 0.5AU and r, 



out 



gem 



1.5AU, where r is the radial distance 



from the central star. This disk model is the standard disk mode 



and 50% more massive than the minimum-mass disk (jHayashi 



for s olar system formation 



198ll ). In the oligarchic 



growth model, the mass of a protoplanet M is given by the isolation mass 

3/2 



M iso ~ 2vra6S = 0.16 



'Si 



lAu) M «' 



(2) 
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where a is the semimajor axis, b is the orbital separation between adjacent protoplanets, 
b = 6/r H , r H is the Hill radius r H = (2M isr ,/ 3M^) 1 / 3 g of the proto planet, M is the mass 



of the central star, and M® is Earth mass (IKokubo fc Ida 



20001). We set the orbital 



separation of protop 



Kokubo fc Ida 



2000 



anets as b = 10 that is the typical value in iV-body simulations (e.g., 



20021 ). The initial eccentricities e and inclinations i of protoplanets 



are given by the Rayleigh distribu tion with dispersions (e 2 ) 1//2 = 2(i 2 ) 1 / 2 = 0.01 (the unit of 



i is radian) ( llda fc Makino 



1992J). We set the bulk density of protoplanets as p = 3 gem 



The initial protoplanet system has the number of protoplanets n = 16 in r- m < a < r out , 
total mass M tot ~ 2.3M e , specific angular momentum j ~ 0.95j e , and mean semimajor 
axis a = j 2 /GM & = 0.91AU, where j® = i/GM « e and a® is the semimajor axis of Earth. 
For each accretion model, we perform 50 runs with different initial angular distributions of 
protoplanets. 



2.2. Orbital Integration 



The orbits of protoplanets are calculated by numerically integrating the equations 
of motion of protoplanets. We set the mass of the central star equal to solar mass. For 
numerical integration , we use the modified Hermite scheme for planetary iV-body simulatio n 



(IKokubo et al. 



1991 



1998 



Kokubo &: Makino 



20041 ) with the hierarchical timestep ( IMakino 



). For the c alculation of mutual 



GRAPE scheme (INitadori et al. 



gravity among protoplanets, we use the Phantom 



20061 1 . The simulations follow the evolution of protoplanet 



systems for 3 x 10 s years until only a few planets remain. 
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2.3. Accretion Condition 



During orbital integration, when two protoplanets contact, a collision occurs. We 
define an impactor/target as a smaller /larger one of two colliding bodies. We obtained 



a realistic accretion conditio n of protop 



The standard SP H method (IMonaghan 



equation of state ( jMelosh 



anets 



1992 



jy performing SPH collision simulations. 



Canup 



2004 ) was used with the Tillotson 



19891 ) . We assumed differentiated protoplanets with 30% core 



and 70% mantle in mass. Protoplanets were represented by 20,000 particles in most runs 
and 60,000 particles in high- resolution runs. We systematically varied the mass ratio of the 
impactor and target as M im p/Mt ar = 1, 2/3, 1/2, 1/3, 1/4, 1/6, and 1/9, their total mass 



as M co i = M imp + M tar = 0.2-2M e , and the impact velocity and angle as v 



imp 



1.0-3.0U, 



at 0.02 or 0.2v esc intervals and 9 = 0°-75° at 15° intervals, where v esc is the mutual surface 
escape velocity v esc = [2GM col / (i?j mp + i? tar )] 1//2 , -Ri mp and i? tar are the radii of the impactor 
and target, and 9 is the angle between the surface normal and the impact trajectory {9 = 0° 
for a head-on collision and 9 = 90° for a grazing encounter). Based on the results of SPH 
collision simulations, we derived an empirical formula for the critical impact velocity, below 
which a collision leads to accretion, in terms of the masses of the impactor and target and 
the impact angle as 



Cl 



Mi 



tar 



Mi 



imp 



[1 -sin#) 5 / 2 + 



c 2 



Mi 



tar 



Mi 



imp 



+ c 3 (l-sinfl) 5 / 2 + C4, (3) 



M col J V M col 

where c\ = 2.43, c 2 = —0.0408, c 3 = 1.86, and c 4 = 1.08 are numerical constants. The 
critical impact velocity is a decrease function of sin 9 as seen in Figure [TJ For collisions 
of equal-mass protoplanets, this conditio n agrees well with the results of SPH collision 



simulations by 



Agnor fc Asphaug 



( 120041 ) . The reason for this agreement is that we 
performed almost the same collision simulations by essentially the same method, which 
confirms the robustness of the results. The details of the collision simulations and the 



accretion condition will be presented in a separate paper. In the realistic accretion model 
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we use this accretion condition. In accretion, the position and velocity of the center of mass 
and the angular momentum are conserved. 

If the impact velocity v = \v\ = \/v 2 + v\ is higher than v cr , two colliding bodies are 
bounced with the rebound velocity v', v' n = and v[ = max(t> t , f esc ), where subscripts n and 
t mean normal and tangential components of the velocity, respectively. The change in spin 
angular momentum brought about by hit-and-run collisions is not taken into account since 
it is usually much smaller than that caused by accretionary collisions. These prescriptions 
follow the results of the collision experiments. 

We assume that, initially, protoplanets have no spin angular momenta. We track the 
spin angular momentum of planets resulting from accretion. When accretion occurs, the 
orbital angular momenta of two colliding bodies about their center of mass and their spin 
angular momenta are added to the spin angular momentum of a merged body. The spin 
angular velocity of planets is computed by assuming a sphere of uniform density, in other 
words, the moment of inertia / = (2/5) MR 2 , where M and R are the planetary mass and 
radius, respectively. 



3. RESULTS 



W e compare the results of the realistic accretio n model to those of the perfect accretion 



model (IKokubo et al. 



2006 



Kokubo fc Ida 



20071 ). The statistics of collision parameters 



and basic dynamical properties of planets derived from 50 runs for each accretion model are 
summarized in Tables [T] and [2j 
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3.1. Statistics of Collision Parameters 

We record the collision parameters in all runs. The numbers of hit-and-run and 
accretionary n acc collisions are summarized in Table [TJ In the realistic accretion model, the 
total number of collisions in 50 runs is 1211. In a run, the average numbers of hit-and-run 
and accretionary collisions are (nhar) = 11-8 and (n &cc ) = 12.4, respectively, which means 
on average 49% of collisions results in hit-and-run. Note that the average number of 
accretionary collisions for the realistic accretion model is almost the same as that for the 
perfect accretion model. This suggests that the growth timescale of planets is independent 
of the accretion model. We will discuss this later. 



On the giant impact stage, the last one or a few 



spin angular momentum of planets (lAgnor et al. 



1999 



arge collisions determ ine the final 



Kokubo fc Ida 



20071 ). In this case, 



the spin angular velocity is estimated as 

u^^g^Y'ismeY 2 ^ 2 ) 1 ^ (4) 

where 7 = M imp /M coh g(j) = 7 (1 - 7) [7 1 / 3 + (1 - 7 ) 1/3 ] 1/2 , and v = v/v csc . The critical 
spin angular velocity for rotational instability u cr is given by 

1/2 / \ 1/2 



0J r 



f gm\ 1 ( p y'\ , 



where p is the material density. 

In Figure (TJ we plot the scaled impact velocity v against sin 9 for all collisions in all runs 
for the realistic accretion model. We show that only collisions with small v and/or sin# can 
lead to accretion due to the realistic accretion condition. For collisions forming Earth-sized 
planets (M > M©/2), the RMS values of the collision parameters are {g{^) 2 ) 1 ^ 2 = 0.24, 
(sin 2 #) 1/2 = 0.62, and (z/ 2 ) 1 / 2 = 1.19 for the realistic accretion model, and (g^) 2 ) 1 ^ 2 = 0.25, 
(sin 2 ^) 1 / 2 = 0.71, and (z/ 2 ) 1 / 2 = 1.40 for the perfect accretion model. The RMS impact 
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parameter (sin 2 ^) 1 / 2 for the realistic accretion model is smaller than that for the perfect 
accretion model since small impact parameters are selective in the realistic accretion model 
Also, for the same reason, the RMS impact velocity v for the realistic accretion model is 
slightly smaller than that for the perfect accretion model. Using these RMS values, we 
obtain the typical spin angular velocity resulting from an accretionary collision as 2.0 hr _1 
and 2.8 hr _1 for the realistic and perfect accretion models, respectively. 



3.2. Statistics of Basic Dynamical Properties 

The average values with standard deviations for basic dynamical properties of planets: 
the number of planets, n, the number of Earth-sized planets with M > M®/2, um, the 
number of planets in r- m < a < r out , n a , the in-situ accretion efficiency f a = M a /M tot where 
M a is the total mass of planets in r in < a < r out , and the growth timescale T grow defined as 
duration of accretion for each planet that undergoes at least two accretionary collisions are 
summarized in Table [TJ The mass and orbital elements (semimajor axis a, eccentricity e, 
and inclination i) of the largest and second-largest planets are summarized in Table [TJ 

First, we find no substantial differences between the results of the two accretion models. 
In both models, a typical resultant system consists of two Earth-si zed planets and one 



2006|). 



or two smaller planets that are as large as the initial protoplanets ( IKokubo et al. 
The two Earth-sized planets tend to form inside the initial distribution of protoplanets. 
The resultant system for the realistic accretion model tends to be slightly wider than 
that for the perfect accretion model, which can be seen as having slightly larger (n) and 
smaller (n a ) and (/ a ). This is because the realistic accretion model experiences more close 
encounters between protoplanets that diffuse the system. It should be noted that though 
about half of collisions do not lead to accretion, this accretion inefficiency does not lengthen 
the growth timescale by a factor of two or more, relative to the perfect accretion model as 
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expected by lAgnor fc Asphaug I (120041 ) . The growth timescale for the realistic accretion 
model is only slightly longer than that for the perfect accretion model. This is because 
even though collisions do not lead to accretion, the colliding bodies stay on the colliding 
orbits after the collision and thus the system is unstable and the next collision takes place 
shortly. In fact, the mean accretionary collision time for the realistic accretion model is 
T co i = 11.1 ± 29.3 x 10 6 years, which is as long as the mean collision time for the perfect 
accretion model T co \ = 10.3 ± 30.3 x 10 6 years. On the other hand, in the realistic accretion 
model, the mean collision time after a hit-and-run collision is short as T co \ = 3.6 ±9.1 x 10 6 
years, where 38% of collisions is for the same pair and its mean collision time is even shorter 
as T col - 2.1 ± 5.6 x 10 6 years. 

Figure [2] shows the average masses of the largest and second-largest planets against 
their average semimajor axes. We find that there are no differences in either mass or 
semimajor axis of the planets for the realistic and perfect accretion models. The largest 
planets with M ~ 1.2M e tend to form around (a\) ~ 0.8 AU, while the second-largest 
planets with M ~ 0.7M® is widely scattered in the initial protoplanet region. We also 
find no difference in their eccentricities and inclinations, and those are ~ 0.1. These 
eccentricities and inclinations are an order of magnitude larger than the proper eccentricities 
and inclinations of the present terrestrial planets. Some damping mechanism such a s 



gravitational drag (dynamical fr iction) from a dissip ating gas disk (jKominami fc Ida 



or a residual planetesimal disk (lAgnor et al. 
stage. 



2002) 



19991 ) is necessary after the giant impact 



3.3. Statistics of Spin 



In 50 runs of the realistic and perfect accretion models, we have 128 and 124 planets 
that experience at least one accretionary collision, respectively. The average values of each 
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component of the spin angular velocity oj of the planets and its dispersion o are summarized 
in Table |2] together with the Root-Mean-Square (RMS) spin angular velocity (uj 2 ) 1 / 2 and 
the spin anisotropy parameter (3 = (oof)/ (oj 2 ). 

The spin angular velocity averaged in mass bins against mass is shown in Figure [3k. 
We show clearly that the average angular velocity is almost independent of mass for both 
accretion models. For the perfect accretion model, the average values are as high as the 
critical angular velocity u cr . The se are natural outcome s for the giant impact stage under the 



assumption of perfect accretion ( iKokubo fc Ida 



20071 ). The RMS spin angular velocity for 



the realistic accretion model is about 30% smaller than that for the perfect accretion model. 
This is because in the realistic accretion model, grazing and high-velocity collisions that 
have high angular momentum result in a hit-and-run, while nearly head-on or low-velocity 
collisions that have small angular momentum lead to accretion. In other words, small 
angular momentum collisions are selective in accretion. Thus, the accretion inefficiency 
leads to slower planetary spin, compared with the perfect accretion model. Indeed, the 
RMS spin angular velocity is almost consistent with the estimation by equation (j4j) in 
section 13.11 The RMS spin angular velocity slightly larger than the estimation based on 
a single dominant accretionary collision is due to the contribution of other non-dominant 
accretionary collisions. 

We confirm that each component of oj follows a Gaussian distribution. We perform a 
Kolmogorov-Smirnov (K-S) test to confirm the agreement with a Gaussian distribution. For 
the distribution of each component, we obtain sufficiently high values of the K-S probability 
Qks > 0.5 in both accretion models. 

In Figure |3b, we show the obliquity distribution with an isotropic distribution, 



nde = (1/2) sinede. We find that the obliquity ranges from 0° t o 180° and fo 



isotropic distribution nde = (1/2) sinede, which is consistent with 



Agnor et al. 



lows an 



(119991 ) and 
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Kokubo fc Ida I (120071 ). By a K-S test, we obtain high K-S probabilities of 0.3 and 0.9 
for the realistic and perfect accretion models, respectively. This is also confirmed by the 
spin anisotropy parameter /3 ~ 1/3 in Table |2J The isotropic distribution of e is a natural 



outcome of giant impacts since the impacts are three-dimensional and add equal. 



contributions to each component of the spin angular momentum ( jKokubo fc Ida 



y random 



2007J). 



For an Earth mass planet, the RMS spin angular velocity 2.35 hr _1 of the realistic 
accretion model corresponds to the spin angular momentum of 9.5 x 10 41 g cm 2 s -1 , which 
is 2.7 times larger than the angular momentum of the Earth-Moon system. So the angular 
momentum of the Earth-Moon system is not a typical value of the realistic accretion model 
but it is reasonably within the distribution of u>. 



4. Summary and Discussion 

We have investigated the basic dynamical properties of the terrestrial planets assembled 
by giant impacts of protoplanets by using iV-body simulations. For the first time, we 
adopted the realistic accretion condition of protoplanets obtained by the SPH collision 
experiments. The basic dynamical properties have been studied statistically with numbers 
of TV-body simulations. For the standard protoplanet system, the statistical properties of 
the planets obtained are the following: 

• About half of collisions in the realistic accretion model do not lead to accretion. 
However, this accretion inefficiency barely lengthens the growth timescale of planets. 

• The numbers of planets are (n) ~ 3—4 and (um) — (n a ) ~ 2. The growth timescale 
is about 6—7 x 10 7 years. The masses of the largest and second-largest planets 
are (Mi) ~ 1.2M® and (M2) ~ 0.7M®. The largest planets tend to form around 
(a^ ~ 0.8AU, while a 2 is widely scattered in the initial protoplanet region. Their 
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eccentricities and inclinations are ~ 0.1. These results are independent of the 
accretion model. 

• The RMS spin angular velocity for the realistic accretion model is about 30% smaller 
than that for the perfect accretion model that is as large as the critical spin angular 
velocity for rotational instability. The spin angular velocity and obliquity of planets 
obey Gaussian and isotropic distributions, respectively, independently of the accretion 
model. 

We confirm that except for the magnitude of the spin angular velocity, the realistic 
accretion model gives the same results as the perfect accretion model. This agreement 
justifies the use of the perfect accretion model to investigate the basic dynamical properties 
of planets except for the magnitude of the spin angular velocity. 

In the present realistic accretion condition, we do not consider the effect of the spin on 
the accretion condition that would potentially change collisional dynamics. It is difficult 
to derive an accretion condition in terms of spin parameters by SPH collision simulations 
since the collision parameter space becomes huge and it is almost impossible to cover all 
parameter space. The fragmentation of planets is not taken into account, either. Including 
the fragmentation may be able to further reduce the spin angular velocity of planets by 
producing unbound collisional fragments with high angular momentum. Furthermore, the 
collisional fragments can potentially alter orbital dynamics of planets through dynamical 
friction if their mass is large enough. The high eccentricities and inclinations of planets 
may be damped by dynamical friction from the collisional fragments. In order to take into 
account these effects and make the model of terrestrial planet formation more realistic, 
we plan to perform iV-body simulations for orbital dynamics and SPH simulations for 
collisional dynamics simultaneously in a consistent way in future work. 
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sin0 

Fig. 1— Scaled impact velocity v is plotted against the scaled impact parameter sin 6 
for accretionary (filled circles) and hit-and-run (open circles) collisions in 50 runs of the 
realistic accretion model. The solid curves show the scaled critical impact velocity v CT /v esc 
for M imp /M tar = 1/10, 1/2, and I. 



Table 1. Basic Dynamical Properties of Final Planets 
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Table 2. Spin Parameters of Final Planets 



accretion model (wx)(br _1 ) (w w )(hr _1 ) (^(hr -1 ) ^(hr" 1 ) ^(hr" 1 ) ^(hr" 1 ) {uj 2 ) 1 ' 2 ^ 1 ) f3 



realistic 0.05 0.08 0.07 1.32 1.26 1.49 2.35 0.40 

perfect 0.42 -0.09 -0.16 2.15 2.17 2.03 3.69 0.31 
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Fig. 2. — Average semimajor axes and masses of the largest (filled symbols) and second- 
largest (open symbols) planets for realistic (circle) and perfect (square) accretion models. 
The error bars indicate 1-er. 
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Fig. 3. — Left: Average spin angular velocity of all planets formed in the 50 runs of the 
realistic (circle) and perfect (square) accretion models is plotted against their mass M with 
mass bin of 0.1M e . The error bars indicate l-a and the dotted line shows u cr . Right: 
Normalized cumulative distributions of e for the realistic (solid curve) and perfect (dashed 
curve) accretion models with an isotropic distribution (dotted curve). 



